Criticality, the area law, and the computational power of PEPS 
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The projected entangled pair state (PEPS) representation of quantum states on two-dimensional 
lattices induces an entanglement based hierarchy in state space. We show that the lowest levels 
of this hierarchy exhibit an enormously rich structure including states with critical and topological 
properties as well as resonating valence bond states. We prove, in particular, that coherent versions 
of thermal states of any local 2D classical spin model correspond to such PEPS, which are in 
turn ground states of local 2D quantum Hamiltonians. This correspondence maps thermal onto 
quantum fluctuations, and it allows us to analytically construct critical quantum models exhibiting 
a strict area law scaling of the entanglement entropy in the face of power law decaying correlations. 
Moreover, it enables us to show that there exist PEPS within the same class as the cluster state, 
which can serve as computational resources for the solution of NP-hard problems. 

PACS numbers: 



The concept of entanglement plays a central role in 
both fields of quantum information theory and of strongly 
correlated systems. In quantum information theory it 
lies at the heart of many applications and it is viewed 
as a resource for various information processing tasks. In 
condensed matter theory entanglement is one of the roots 
for the notorious complexity of quantum many-body sys- 
tems: its presence necessitates a description within an 
exponentially growing Hilbert space and it is intimately 
connected with many of the fascinating properties which 
quantum matter can exhibit at small temperatures. 

Many fundamental questions arise at the crossing of 
these fields: how is entanglement related to the power 
of quantum computation on the one hand, and the dif- 
ficulties of classical simulations on the other? What is 
the scaling of the entanglement entropy in spin systems, 
its relation to criticality, and the appearance of topo- 
logical quantum order? All these questions can be ad- 
dressed very easily within the framework of so-called 
projected entangled pair states (PEPS) — this is the in- 
tention of this paper. We will see in particular that all 
the above mentioned properties emerge naturally already 
within the simplest classes of PEPS which include cluster, 
toric code and resonating valence bond states. This will 
enable us to settle a recent debate about the relation be- 
tween criticality and entropy scaling, and it allows us to 
find computational resources for the solution of NP-hard 
problems. The central tool of the paper is a general cor- 
respondence between thermal states of classical 2D spin 
models and 2D quantum states with a simple PEPS rep- 
resentation. This correspondence substitutes thermal by 
quantum fluctuations while preserving the nature of cor- 
relations, and it thus maps critical classical onto critical 
quantum states. 

We begin by recalling the PEPS formalism, which was 
introduced in the context of numerical renormalization 
group methods for simulating strongly correlated quan- 



tum spin systems 0,0- PEPS can be viewed as general- 
izations of the AKLT valence bond solids to arbitrary 
lattices and dimensions. Consider an arbitrary connected 
graph where each of N vertices corresponds to a quan- 
tum system, a spin, with d degrees of freedom. A PEPS 
\^) e C is then constructed by (i) assigning to each 
vertex as many virtual spins of dimension D as there 
are adjacent edges, (ii) putting a maximally entangled 
state |/) = Y^iLi N) onto each edge, and (iii) mapping 
the virtual onto the physical spins by applying a map 
P : «) ... (g) ^ C'' at each vertex. Naturally, the 
graph is chosen according to the physical symmetry, and 
although most of the following holds in general we will 
consider square lattices throughout. 

The power of the PEPS formalism is based on two 
points. First, every state has a PEPS representation 
r^]. Hence, with increasing D this representation induces 
a hierarchy in the space of states, from product states 
{D = 1) to more and more entangled ones. Second, 
it appears that many states arising in physics are very 
well approximated by the lower levels of this hierarchy 

. This makes them a powerful variational class for nu- 
merical renormalization group methods on the one hand 
P, Q , and an interesting testbed for all kinds of quantum 
many-body questions on the other . 

Quantum- classical correspondence: Consider a 
classical two-body spin Hamiltonian of the form 
H{ai, ■ ■ ■ ,aN) = J2{i,]) ^(o'l-o-j) with ai ^ 1, . . . ,d and 
respective partition function Z — J2a^^P[~l^-^i'^)] 
inverse temperature f3. From this a corresponding quan- 
tum state can be constructed by using the Boltzmann 
weights as superposition coefficients such that 



\^H.,t3) 



— y 



,o"n) 



(1) 



We will see that \tpH,p) has the following properties: (i) 
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for diagonal observables it gives rise to the same expecta- 
tion values and correlations as the classical thermal state, 
(ii) it has a simple PEPS representation with D — d, (iii) 
it is the ground state of a local quantum Hamiltonian, 
and (iv) when considering asymptotically large systems 
{N oo) the scaling of the entropy of a block of spins 
obeys a strict area law. Whereas (i) is a direct conse- 
quence of the construction, (iii) and (iv) are implied by 
the PEPS parametrization. In order to see the latter we 
rewrite the state as 
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exp 



f3 
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where |+) = X]s=i I*) ^^^"^ ^ diagonal operator 

acting on sites i,i as hij \ai,aj) — h(ai,aj)\ai,<7j). Fol- 
lowing Eq.(j21l we can think of the state \4'H,f3) as being 
constructed from the product state |-|-, • • • , by ap- 
plying (non-unitary) gates exp[— /3/i/2] to all neighbor- 
ing spins. In fact, we may interpret Eq.(|2J) as a quan- 
tum cellular automaton evolution in imaginary time. As 
explained in a nonlocal gate like exp[—Ph/2] can 
be reexpressed by local operations which act addition- 
ally on an auxiliary maximally entangled state. More 
specifically, we take operators V, V' : C^, each 

acting as P\s,k) = \s){(ps\k). Then we obtain indeed 
exp[— /3/2/i] — {V®T") \I) if we choose the vectors fsi^'s' 
such that '^'l^i{^s\k){(p'g,\k) = h{s,s') which is always 
possible, e.g., by a singular value decomposition. Ap- 
plying these gates to all edges leads then to the desired 
PEPS representation. 

As an example consider the ferromagnetic Ising model 
on a 2D square lattice with 
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±1 . 



In this case we can choose ips = f'^ such that {<fs\k) are 
the matrix elements of the square root of the matrix h. 
Applying all the gates gives then rise to a PEPS with 

P = |0)(^ol(^ol(^ol(^ol + . 

Clearly, the expectation values of Pauli Sz operators in 
\'4'H,i3) equal the classical expectation values. In the 
quantum case, however, we do not only have diagonal 
observables, but also non-diagonal ones like Sx- Surpris- 
ingly, their expectation values are determined by classical 
ones as well, like 



E 
E 



H{<Ti ,CT2,...)+-ff( — (Ti ,(T2,...) 



0'1,(T2, 



where the last term is a local 5-body expectation value 
in the classical Gibbs state. In general, every local ex- 
pectation value in the quantum state corresponds to a 




Figure 1: The entropy of a block of spins in a PEPS scales like 
the perimeter of the block: the Schmidt rank of the reduced 
density operator of the considered block of spins is bounded 
above by the product of the Schmidt ranks of the broken 
bonds. 



local expectation value in the classical state, where the 
region the observable acts on is enlarged at most by the 
interacting neighborhood. 

Before continuing it should be noted that several re- 
sults related to the above classical-quantum correspon- 
dence can be found in the literature: a connection 
between so-called Rokshar-Kivelson points and classi- 
cal stochastic models was recently made in Q and 
between Hamiltonians and rapidly mixing reversible 
Markov chains in 0. In a generalization of the 
AKLT-state on 2D lattices was considered and demon- 
strated that it can be mapped onto a classical vertex 
model. The PEPS formalism, discussed in the present 
paper, provides a very natural framework for describing 
and generalizing those results. 

Criticality and the area law: Recently a lot of atten- 
tion has been devoted to the scaling of the entanglement 
entropy 0, O [13, IH Ell- That is, given a ground 
state, how does the entropy of a contiguous subsystem 
scale with the size of the latter? Originally appearing 
in the context of black holes the renewed interest in this 
question comes from the investigation of quantum phase 
transitions and the quest for powerful ansatz-states for 
the classical simulation of quantum systems. In ID it 
is known that critical states corresponding to a confor- 
mal field theory exhibit a logarithmic divergence of the 
entropy S ^ logL {L being the length of the subsys- 
tem), whereas there seems to be a saturation for all 
non-critical systems llll. In I? > 1 dimensions quasi- 
free systems of Bosons [l^ and Fermions [ISj have been 
studied. Whereas in the non-critical (gapped) bosonic 
case there is a strict area law S = 0{L^~^) (for a cube 
with edge-length L), this is violated in the gapless case 
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of Fermions, where S ~ L^~-^ logi. This naturally rises 
the question about a one-to-one correspondence between 
criticality and a violation of the area law. The PEPS for- 
malism together with the above classical-quantum corre- 
spondence enables us now to answer this question in the 
negative in a very simple way: consider the classical 2D- 
Ising system, which is known to become critical in the 
thermodynamic limit at /3c = ^(l + \/2)- The corre- 
sponding quantum state will then have exactly the 
same correlations {S'^SD ~ l/v^, which now reflect crit- 
ical quantum rather than thermal fluctuations. In spite 
of this, the state obeys a strict area law bound S < 4L 
due to the PEPS representation with D = 2 since the en- 
tropy is solely generated by breaking 4L entangled bonds 
at the boundary (see figure 1). In the following we will 
show that, unlike the cases in |Q, the power law de- 
cay of correlations is not a consequence of long range 
interactions. Rather {ipfjj is the ground state of a local 
Hamiltonian. 

Parent Hamiltonians: Every PEPS with finite D is 
(on a sufficiently large lattice) the ground state of a lo- 
cal Hamiltonian. The standard construction of such par- 
ent Hamiltonians identifies projectors onto null spaces of 
reduced density operators with the interaction terms in 
the Hamiltonian Here we will follow a different 

approach, related to the one described in 0, 0|, which 
is adapted to PEPS constructed from classical models 
and allows us to prove uniqueness and existence of a gap 
above the ground state energy for P < (3c- 

Consider any ergodic local Markov process obeying de- 
tailed balance and converging to the equilibrium distribu- 
tion of the classical model at inverse temperature (3 (e.g., 
by using Metropolis Monte Carlo or general Glauber dy- 
namics ^3). The stochastic transition matrix corre- 
sponding to the Markov process can be written as a sum 
M(/3) = J2k^k{l3) where each Mfc(/3) acts locally and 
obeys detailed balance. The latter requires 

p,(/3)[Mfe(/3)]^^^=Pb(/3)[Mfc(/3)]^^^ 

with PaiP) — e'xjp{—f3H{a))/Z and a,b denoting a par- 
ticular configuration of the N spins. In matrix notation, 
this is equivalent to imposing that all the matrices 

Pfc(/3) = e"^ ^ w) '^'^ MkiP) ^ w) 

are symmetric. Obviously, the operator Pk{P) is sym- 
metric and has exactly the same eigenvalues as J2k ^kiP) 
since they are connected by a similarity transformation. 
Furthermore, all Pk{(3) are local operators if the Mk{(3) 
were local. Note that 1 — J2k^^k{(3) only has non- 
negative eigenvalues with the equilibrium distribution 
corresponding to eigenvalue 0. We thus define the Hamil- 
tonian H{f3) = 1 — J2k Pk{P) > which is a sum of local 
operators such that by construction jV'/s) is the ground 
state of H{P). Moreover, H{f3) is gapped iff the stochas- 
tic matrix M(/3) has a gap. In this way the gap in the 



quantum Hamiltonian corresponds to the rate of conver- 
gence to equilibrium of the Markov process. In fact, the 
existence of a gap in M(/3) for (3 < (3c was proven in 0| 
for a class of models including the 2D Ising model. At 
precisely the critical point, Monte Carlo methods exhibit 
a slowing down, leading to a gapless critical quantum 
Hamiltonian. In fact, power law decaying correlations 
imply that the Hamiltonian has to be gapless [l^ . 

Computational power of PEPS: In this section we will 
treat the PEPS as a resource for computational tasks. 
Given a source which produces a specific state in an effi- 
cient manner, together with the ability of performing ar- 
bitrary local measurements, what kind of computational 
problems can we solve efficiently? This question is clearly 
inspired by the cluster state computational model |l8j | . In 
fact, it was shown in that the cluster state is a PEPS 
with D = 2. Moreover, it is known to be a resource state 
for universal quantum computation, i.e, it enables us for 
instance to solve a typical NP problem — factorization — 
by merely performing local measurements. Since PEPS 
with D — 1 are product states, _D = 2 is in fact the 
simplest class in which useful resources can be expected. 
Exploiting the above formalism it is now simple to show 
that there are other powerful resource states within this 
class, which even enable the efficient solution of NP-hard 
problems. In order to see this, note that given a quantum 
state IV'/?) which corresponds to a classical Hamiltonian 
H{a) we can efficiently determine expectation values in 
the classical Gibbs state by performing local (diagonal) 
measurements on {ipp}- For /? — > oo we can for instance 
measure the ground state energy This is in partic- 
ular true for the D = 2 PEPS corresponding to a two- 
dimensional Ising spin glass within a magnetic field. This 
task was, however, shown to be an NP-hard problem in 
[20| together with the case of 3D spin glasses without 
magnetic field. Similarly, the determination of the parti- 
tion function of the Potts model (£> > 2) is known to be 
^P-hard and tightly connected to hard problems in knot 
theory. As the task of calculating expectation values of 
PEPS can be done by contracting a network of tensors 
arranged on a square lattice 0, the above arguments 
prove that such a contraction of tensors is in general a 
NP-hard problem in the number of tensors jl^l • 

The approach of encoding the solution to an NP-hard 
problem into a quantum state is reminiscent of adiabatic 
quantum computing which, however, deals typically 
with ground states of ID albeit non-local Hamiltonians. 
In fact, one possible way of generating PEPS would be by 
adiabatic means with the usual caveat concerning the gap 
of the system. However, as in the case of the cluster state, 
there might be better ways of generating these states 
since after all we have an efficient local parametrization — 
the above observation makes the generation of PEPS a 
highly interesting problem. 

Let us finally show that two other classes of states, im- 
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portant for quantum information and condensed matter 
theory, are contained within small- Z? PEPS as well: 

Toric code states introduced in the context of quan- 
tum error correction are very interesting as they exhibit 
nontrivial topological properties [i^. In the case of an 
infinite square lattice, the toric code is the ground state of 
a Hamiltonian consisting of local commuting projectors, 
each of them annihilating the ground state. The state 
can again be written in terms of the (zero-temperature) 
Boltzmann weights of a classical statistical model 



lim exp 




where denotes the i'th plaquette in the lattice with 
the spins on the edges. We can again represent the com- 
muting nonlocal gates by introducing entangled auxiliary 
degrees of freedom and applying local operations. More 
specifically, the 4-qubit gate can be implemented by dis- 
tributing the states |/) between the 4 qubits, followed 
by local projections of the form |0)(5'+| + where 
|vl/±) = |00)±|ll). 

Implementing this map on all plaquettes of the lattice, 
we obtain projectors of the form 



Pe = |0)(*+|l2(*+|34 
Po = |0)(«'+|i4(*+|23 



|l>(*"|l2(*-|34 
|l)(*-|l4(*-|23 



where Pg, Po act on the even and odd sites of the bipar- 
tite lattice respectively (the labels 1..4 denote the virtual 
qubits in clockwise order). Hence it is again a simple 
PEPS with D = 2 exhibiting nontrivial topological be- 
havior. In the case of PEPS with finite Z?, it is indeed 
always simple to calculate the topological entropy defined 
in [2^1 explicitly — the PEPS formalism seems to provide 
a promising avenue for generating other states exhibiting 
those fascinating properties. 

Resonating valence bond states (RVB) have been stud- 
ied extensively in the context of strongly correlated sys- 
tems [i^. These states exhibit topological quantum or- 
der and do not seem to have any classical statistical 
model associated to them because the wave function con- 
tains negative weights. For the case of simplicity, let 
us consider the simplest RVB state which is the equal 
weight superposition of all possible coverings of singlets 
over nearest neighbors on a square lattice. It can eas- 
ily be checked that this RVB is equivalent to the PEPS 
defined by 

P = |0> ((0222| + (2022| + (2202| + (2220|) 
+ |1) ((1222| + (2122| + (2212| + (2221|) 

acting on virtual singlets of the form \S) = |01) — |10) -|- 
|22) distributed between all nearest neighbors. Interest- 
ingly, we need 13 = 3 in this case, and again the area 



law is automatically proven (i.e., the entropy of a block 
of spins scales like the boundary) . In a similar way, RVB 
with singlets distributed beyond nearest neighbors can 
easily be constructed. 

In summary we found that already the lowest lev- 
els of the PEPS hierarchy exhibit an enormously rich 
structure — they contain highly interesting states for 
quantum information (e.g., cluster and toric code states) 
as well as for condensed matter theory (e.g., critical and 
RVB states). This makes them an interesting variational 
class for numerical methods and a rich testbed for quan- 
tum many-body questions. Based on a classical-quantum 
correspondence we were able to find critical quantum 
models whose entropy scaling contrasts with the one for 
Fermions (and corresponding spin models jl^). More- 
over, it yielded a local description of simple PEPS en- 
coding the solution of NP-hard problems. 

Work supported by the Gordon and Betty Moore 
Foundation, European projects, der Bayerischen Staat- 
sregierung and MEG MTM-2005. 
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